A Wigner Surmise for Hermitian and Non-Hermitian Chiral Random Matrices 
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We use the idea of a Wigner surmise to compute approximate distributions of the first eigenvalue 
in chiral Random Matrix Theory, for both real and complex eigenvalues. Testing against known 
results for zero and maximal non-Hermiticity in the microscopic large- iV limit we find an excellent 
agreement, valid for a small number of exact zero-eigenvalues. New compact expressions are derived 
for real eigenvalues in the orthogonal and symplectic classes, and at intermediate non-Hermiticity 
for the unitary and symplectic classes. Such individual Dirac eigenvalue distributions are a useful 
tool in Lattice Gauge Theory and we illustrate this by showing that our new results can describe 
data from two-colour QCD simulations with chemical potential in the symplectic class. 
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1. Motivation. Probably one of the most used predic- 
tions of Random Matrix Theory (RMT) is the so-called 
Wigner surmise (WS) describing the universal repulsion 
of energy levels in many systems in nature, including neu- 
tron scattering, quantum billiards and elastomechanical 
modes in crystals For large matrices, the nearest- 
neighbour (nn) spacing distribution pW(s) is universal 
and only depends on the repulsion strength which takes 
discrete values (3 = 1,2,4 for the three classical Wigner- 
Dyson (WD) ensembles. It can be computed with sur- 
prising accuracy using 2x2 matrices, which is the WS. 
Although simple arguments discussed in [2J lead to this 
rule for (3—1, such an approximation is by no means 
obvious. 

The extension from WD to non- Hermitian RMT intro- 
duced long ago by Ginibre [3| has become a very active 
field in the past decade, in particular due to applications 
in open quantum systems, see Q for references and other 
applications. Here the spacing is known only for the class 
with broken time-reversal (J3 = 2) and has been applied 
in Lattice Gauge Theory (LGT) ||. However, a simple 
surmise based on 2 x 2 matrices does not work here. 

In this paper we investigate the existence of a surmise 
for the smallest eigenvalue in chiral RMT and its non- 
Hermitian extensions. These have become relevant due 
to applications in Quantum Chromodynamics (QCD) ini- 
tiated by [6] and extended to non-Hermitian QCD at 
finite quark chemical potential fi 0- QCD at strong 
coupling is a notoriously difficult theory, and the chiral 
RMT approach has become an important tool for LGT 
with exact chiral fermions For non-Hermitian QCD 

the complex action hampers a straightforward LGT ap- 
proach, see [loj for a recent discussion and references. 
Here RMT predictions remain possible for various quan- 
tities 0,03, El. 

In this paper we will show that an excellent approxi- 
mation for the 1st non-zero eigenvalue is possible using 
a simple 2 x (2 + v) matrix calculation, capturing the 



repulsion of a small number v of zero eigenvalues. Being 
localised and non-oscillatory the 1st eigenvalue is much 
more suitable for LGT than the spectral density, compare 
e.g. [§] and Our surmise fills some gaps in predic- 
tions for real eigenvalues in the orthogonal and symplec- 
tic classes ((3 — 1,4) [ljj . where until very recently nu- 
merically generated RMT had to be used for comparison 
[lf3 |. We also provide new predictions for intermediate 
non-Hermiticity and test them against QCD-like LGT 
data from [17] . This further completes the picture, com- 
pared to previous approximations [14] ((3 — 2) based on a 
Fredholm determinant expansion [18j |. and exact results 
at maximal non-Hermiticity [lj| ((3 — 2,4). 

2. Level spacing in the WD class. We recall here the 
success of a WS for Hermitian, and its failure for non- 
Hermitian, WD ensembles. The WD partition function 
for an TV x TV Hermitian matrix H with real, complex or 
quaternion real entries is given terms of eigenvalues by 
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The Jacobians of the corresponding ensembles which are 
called GOE, GUE and GSE {(3= 1,2 and 4) include the 
Vandermonde determinant, Ajv(A) = Yik>i(^k — A/). 

The large-TV nn spacing in the bulk of the spectrum 
can be computed approximately from TV = 2 (WS) by 
inserting S(\\i — Aa| — s) in Z- 



WD' 



Pwl( S ) = a ex PM/3 • 
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The constants ap, bp follow from fixing the norm and first 
moment to unity (see e.g. in 1]). The latter can always 
be achieved from dssp(s) — m by rescalingp^^(s) = 
mp^(ms). This fixes the scale compared with TV = oo. 

The exact result p™(s) is cumbersome, given in terms 
of an infinite product of eigenvalues of spheroidal func- 
tions (e.g. in [2]), the 5th Painleve transcendent Q, or 
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FIG. 1: Left: surmise p™g(s) (red) vs exact result (dashed 
blue) in 0- Ri S ht: PaL( r ) for ^ = 2,3,4,20 (red to blue). 

combining a Taylor series with coefficients given by sums 
over permutations and Dyson's asymptotic expansion in 
a Padc approximation [201 ] . This is compared to the sur- 
mise Eq. ((21) in Fig. Q]lcft. In Table U we give the root 
of the integrated square deviation for later comparison, 
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The non-Hermitian WD ensembles are defined by drop- 
ping the Hermiticity constraint in Eq. ([T]) left [3j. We 
only display the complex eigenvalue representation for 
(3 = 2 and 4 and their Jacobians Jp{z) computed in 
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J 2 {z) = \A N (z)\ 2 , Mz) = A 2N {z,z*)]J(z j -z*). 

3=1 

For (3 — 2 the spacing is obtained from an N = 2 surmise 
by inserting 5{\zi — z%\ — r) in Zq!, and putting one 
eigenvalue at the origin. The exact spacing for any N 
obtained in 1211 uses translational invariance in the bulk 
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In Fig. [T] right we compare N — 2 with increasing N, 
all curves having norm and first moment 1. Clearly a 
surmise does not work for the (3 = 2 Ginibre ensemble 
(6 « 0.18), as previously noted in [2l|. For (3 = 4 and 1 
the spacing is currently unknown. 

3. First eigenvalue in chiral RMT. The chiral ensem- 
bles with real eigenvalues called chGOE, chGUE, and 
chGSE are defined in terms of N x [N + v) rectangular 
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TABLE I: Deviation Eq. Q in units 10 3 between approxi- 
mate N = 2 and exact large- iV results (<5ws from |20|). 



matrices W with real, complex or quaternion real ele- 
ments without further symmetry restrictions. Switching 
to positive eigenvalues A^ > of the Hermitian Wishart 
(or covariance) matrix WW^ we obtain 
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Here Nf massless flavours can be added by shifting d — > 
d + Nf. The gap probability E^(s) that the interval 
(0, s) is void follows by integrating in Eq. © from s to 
oo. For N = 2 we obtain 

y* oo 

EW(s)~ dxdy[(x + s){y + x + s)} d c-^ s+x '>-yy' 3 7 (7) 



after shifting variables. The nested integrals can easily be 
evaluated. Note that d = for (3 = 2, v = 0, and /3=l = f. 
These gap probabilities can be computed exactly for any 
TV, and our surmise gives the exact result after rescaling. 

To compare with Dirac operator eigenvalues we have to 
switch variables A 
and thus to s — > s 



3 —r i/j, coming in eigenvalue pairs ±yj, 
1 The distribution of the first positive 



Dirac eigenvalue follows: ptf\s) 

We first list all its known Nf = results in the universal 
microscopic limit for v 6 N in Eqs. ([5]) - (HOD: the chGUE 
for all v [IE SI, the chGOE for v = [3] and odd v 
[H, and the chGSE for v = Q and v > 0. For 
the latter, only a convergent Taylor series is known with 
coefficients a,j(v) given by sums over partitions (see Eq. 
(8) in 24]), much alike for the WS in the WD class, 
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Next, we give examples following our surmise Eq. ([7]) 
where p^(s) is not known in elementary form, filling the 
gaps in Eqs. © - (fTU]) for the first two values of v > 0: 
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The normalisation constants suppressed above easily fol- 
low. However, we cannot set the 1st moment to one as 
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FIG. 2: p^\s) with v = 0, 1, 2 in dashed blue to green for 
(3 = 2 (top left), ,3 = 4 (top right) and (3 = 1 (bottom). The 
N = 2 surmise is in red. Our new result for i/ = 2 at /? = 1 
is compared to a numerical simulation at N = 20 (black dots) . 



in the WD class. The position of pi^ (s) measures the 
repulsion by v exact zero-eigenvalues, containing impor- 
tant information. Thus we fix the N = 2 scale by setting 
the 1st moment equal to the exact one. Without exact 
((3 = 1, even v) or concise ((3 = 4, v > 0) results we in- 
stead fit to the increasing slope of the known microscopic 
density pi^\s), being the first term in the Fredholm ex- 
pansion of the 1st eigenvalue [l8| (see also Eq. (fT8|) ). In 
Fig. [5] we compare approximate to exact 1st eigenvalues 
for small topology v — 0,1,2 and all (3. The deviation 
measured by Eq. ([3]) in Table U increases with v, becom- 
ing visible only for v — 2 (see Fig. [2]). This has to be 
compared to the statistical error in data, see e.g. Fig. [5] 

Note that in chiral RMT the nn spacing also obeys Eq. 
([2]) , but does not follow from an N — 2 surmise [25| ■ 

The non-Hermitian chiral ensembles with /i ^ are 
given in terms of a two-matrix model 11, 2(|. We only 
focus on (3 = 2, 4 here, with their complex eigenvalue 
representations for Nf = reading ll|, l2q 
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The weight w(z) depends on a 
with /i G (0,1]- The limit fi 



1±J£ 
~~2^ 



> b 



> 0, 



leads back to real 



eigenvalues, and at /i = 1 non-Hcrmiticity is maximal. 
The definition of a gap probability on C is not unique 
14 , 1|| . For radial ordering it reads 



EW(r)~H drjrj , 

j=l J r JO 



dBjw{zj)Jp{z) . (14) 



Differentiation yields d r E« 3 \r) = d9 pi P) {re w ), the 
integrated 1st eigenvalue. For (3 = 2 (4) the gap proba- 



bility is given by a Fredholm determinant (Pfaffian) |19| 
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Its matrix elements A^J can be computed recursively for 
any v by differenting the following matrix element [19(: 

(0) _ br 2 h(br 2 )K a (ar 2 ) + ar 2 hibr^K^ar 2 ) 
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This leads to a (3 x (3 determinant (Pfaffian) representa- 
tion for our N = 2 surmise valid for any /i. At /i = 1 

all Fredholm eigenvalues 1 — X^2q n-i are explicitly 
known [19J, providing an exact result for any N as in Eq. 

^ . It contains incomplete Bessel function series (x) 
truncated at power k (= for k < 0) 
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+r 2 (j| fe + fkr 2 )^ +1 (r 2 ) + /^- 1 l(r 2 )^ +2 (r 2 )). 

For (3 = 4 we have the relation = X^k+i with v — > 2v 
[l9j |. In Fig. [3] we compare our surmise to this result, 
truncated at N = 8 because of rapid convergence. Here 
it works better for /3 = 4 than (3 — 2, in contrast to /i = 0. 
Due to angular integration only one scale has to be fixed 
after normalisation, which can be done as in the real case. 




FIG. 3: Integrated 1st eigenvalue at — 1 for i/ = 0, 1,2 (blue 
to green dashes) vs N = 2 (red): (3 = 2 (left) and /3 = 4 (right). 

Next we give a surmise for pi^\re l9 ). In Eq. (fT4| we 
skip the integration over 0\ and differentiate wrt r\. For 
N = 2 we obtain an exact Fredholm expansion 

p^\z)=R[%)- f 1 dtt [ * dtpR^frte*) , (18) 
Jo Jo 
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iy. The 1- and 2-point spectral 



densities are expressed through the kernel of orthogonal 
Laguerre polynomials of norm hj (see [ll| for details) 
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andR^(z,u) = R^(z)R^(u)~\K^(z,u*)\ 2 . For (3 = 
4 we have a Pfaffian of a matrix kernel instead [26j|. An 
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An alternative to Eq. ([18]) is the truncated Fredholm 



FIG. 4: fdO pi 4) (re ie ) (red) vs Lattice data Q3] (blue) with 
volume V — 4 4 , gauge coupling 1.3, /iLat = 0.2 and mass 20 
in Lattice units, using a very large number 10 5 configurations. 



example for p„l (z) is shown in Fig. [5] top right. Here 
two scales have to be fixed: for z we fit to the increase 
of the known microscopic density in the a;-direction, and 
for rescaling 2iV/i 2 = ct 2 to its decrease in the y-direction. 
Since a < 2 for N — 2, we conclude that at large- AT for 
a > 2 pi^ (z) must become symmetric wrt rotation ((3 = 2) 
or reflections wrt the bisector of each quadrant ((3 — 4). 
We have checked this, as well as distributions for < fi < 1 
by generating ensembles of large random matrices. 

4. Lattice data. In two-colour QCD was compared 
to the (3 = 4 microscopic _spectral density in the complex 
plane from chiral RMT [26j. We use the same data here 



but with higher statistics, and refer to [17| for all sim- 
ulation details. Because unimproved staggered fermions 
are used we are in the (3 = 4 class at v = 0. Our Nf — 2 
data are effectively quenched for the smallest eigenvalues 
due to a large mass. In Fig. 0] we compare to the 1st 
integrated eigenvalue, with a — 1.352 being close to max- 



imal non-Hermiticity. No further fits compared to 17 1 
are made. In Fig. [5] we compare LGT data at intermedi- 
ate /^Lat = 0.1 to the angle-dependent surmise Eq. (fl8|) 
by taking cuts. Here the two scales are fitted to the data, 
finding an excellent agreement for a = 0.65. 




expansion in the microscopic large- N limit 18] which was 
successfully applied to the (3 — 2 class [T3]- However, inte- 
grals of higher order terms rapidly become cumbersome. 

5. Conclusions. Conceptually it is possible within chi- 
ral RMT to approximate the 1st eigenvalue distribution 
using a 2 x (2 + v) matrix calculation, for both real and 
complex eigenvalues. It is remarkable that this surmise 
works and captures the repulsion of v zero-eigenvalues. 
We derived new compact expressions for (3 = 1 and 4 
with real eigenvalues for v > 0. Second, we have shown 
that our surmise for (3 — 4 successfully describes SU(2) 
Lattice data, in an intermediate regime for /i ^= where 
no results were previously known. It would be very inter- 
esting to extend our results to the (3—1 non-Hermitian 
chiral class, having both real and complex eigenvalues. 

Support by ENRAGE MRTN-CT-2004-005616 (G.A., 
E.B.), EPSRC grant EP/D031613/1 (G.A., L.S.) and 
DFG grant JA483/22-1 (E.B.) is acknowledged. 



FIG. 5: Top: contour plots for Lattice data as in Fig. g] but 
with /iLat =0.1 (left) vs surmise Eq. (|18[) (right). Bottom: 
cuts through a single peak in x- (left) and y-direction (left). 
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